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THEORY  OF  BOUND  STATES  ASSOCIATED  WITH  n-TYPE 
INVERSION  LAYERS  ON  SILICON* 
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Irvine,  California  92717 


Abstract 

The  ground  state  and  two  excited  states  of  an  electron  bound 
to  a charged  impurity  located  at  the  interface  between  silicon  and 
silicon  dioxide  have  been  investigated.  The  interface  was  taken 
to  be  parallel  to  a (001)  plane  and  a static  electric  field  per- 
pendicular to  the  surface  was  assumed.  Calculations  of  the 
electron  binding  energies  as  functions  of  electric  field  were 
made  for  a bare  charged  impurity  with  screening  by  free  carriers 
neglected.  It  was  found  that  the  binding  energy  increases  with 
electric  field  and  approaches  that  of  the  corresponding  state 
of  a two-dimensional  hydrogen-like  atom  in  the  limit  of  infinite 
electric  field.  The  case  where  the  impurity  is  located  within 
the  oxide  has  also  been  investigated. 
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I . INTRODUCTION 

The  development  of  MOS  devices  has  stimulated  a great 
deal  of  Interest  in  the  basic  properties  of  inversion  layers 
at  semiconductor  boundaries.  A number  of  these  properties  are 
Influenced  by  the  presence  of  impurities  and  their  associated 
electronic  bound  states.  We  have  Investigated  using  a vari- 
ational technique  the  ground  state  and  two  excited  states 
of  an  electron  bound  to  a bare  charged  impurity  located  at 
the  interface  between  silicon  (001)  and  silicon  dioxide.'' 

Theoretical  investigations  of  the  wave  functions  and 
binding  energies  of  donor  states  localized  at  or  near  semi- 
conductor surfaces  without  inversion  layers  have  been  carried 
out  by  Karpushin^^^  and  by  Bell  et  al.^^^  The  effect  of  n-type 
inversion  layers  and  electric  sub-bands  has  been  studied  by  Stern 
and  Howard  for  InAs  and  Sl(OOl)  surfaces.  They  considered 
the  screening  produced  by  the  inversion  layer  charge  carriers, 
but  restricted  their  numerical  results  to  the  two-dimensional 
limit.  In  the  present  paper  we  eliminate  the  latter  restriction, 
but  neglect  the  effect  of  screening. 

For  the  case  of  an  n-type  impurity  at  a semiconductor-oxide 
Interface  where  the  potential  energy  rises  dlscontinuously  to 

infinity,  the  ground  state  is,  in  hydrogen  atom  nomenclature, 

(2) 

the  2p_  level;  ' In  addition,  the  first  excited  states  to 
o 

which  a transition  might  be  observed  in  optical  absorption  are 

^ 4) 

the  3d^j^  and  4d^j^  states.  This  result  is  due  to  Levine'  ' who 
showed  that  a state  cannot  exist  unless  l m is  an  odd  Integer 
and  that  both  t and  m must  differ  by  ±1  in  an  electromagnetic 
transition.  Here  we  consider  only  the  3d^^  excited  states. 


3. 


The  binding  energy  of  the  impurity  ground  state  which  we 
calculate  is  compared  to  the  recent  experimental  data  of 
Hartstein  and  Fowler^®^ 

II.  THEORETICAL  DEVELOPMENT 

We  consider  two  contiguous,  semi-infinite  half-spaces, 
one  of  p-type  silicon  and  the  other  of  Si02)  with  a common 
boundary  parallel  to  a (001)  plane  of  the  silicon.  We 
assume  that  the  potential  energy  of  an  electron  undergoes  a 
discontinuous  Jump  as  the  electron  passes  from  the  silicon 

(S') 

into  the  Si02.  This  discontinuity  in  energy  is  about  3 eV  ; 

however,  we  shall  assume  that  the  discontinuity  is  infinite  and 

that  the  electronic  wave  functions  vanish  at  the  boundary.  In 

the  region  of  the  inversion  layer,  there  is  an  electric  field 

present  which  in  general  is  a complicated  function  of  the  dis- 

( 3 ^ 

tance  from  the  interface.  To  a good  first  approximation'  one 
can  take  the  electric  field  to  be  constant,  a procedure  which  we 
adopt  in  this  paper. 

We  assiuae  that  an  impurity  ion  of  charge  -i-Ze  is  located 
at  the  boundary  between  the  silicon  and  the  Si02.  For  this 
situation,  Karpushin^  has  shown  that  in  the  absence  of  an 
external  electric  field  the  lowest  lying  bound  states  associated 
with  a (001)  Interface  on  silicon  are  primarily  derived  from  the 
energy  ellipsoids  whose  major  axes  are  perpendicular  to  the 
interface.  This  situation  continues  to  hold  in  the  presence 
of  an  electric  field.  We  therefore  consider  only  the  ellipsoids 


just  mentioned.  We  also  neglect  intervalley  interactions  between 
the  two  ellipsoids  perpendicular  to  the  surface.  Finally,  we 
operate  within  the  framework  of  the  effective  mass  approximation. 

In  the  light  of  the  foregoing  discussion,  we  employ  the 
Hamiltonian 


Y + U(r)  + ^ 4-  X ez 

dz 


where  y “ ^ *"  ^ ^ /4Z^m^e  , 

£ is  the  external  electric  field,  m^  and  m^  are  the  transverse 

and  longitudinal  effective  masses,  respectively,  and  and 
are  the  dielectric  constants  of  the  Si02  and  silicon,  respectively. 
In  Eq.  (1),  energy  is  measured  in  units  of  the  effective  Rydberg, 
Ry*  - 2m^e‘*Z^/h2(6j^+€2^^  and  length  in  units  of  the  effective  Bohr 

^ O O 

radius,  a ■ h (€,  +€o)/2Zm  e . For  the  silicon-silicon  dioxide 

O X 2 X 

system,  Ry  - 42.3  meV  and  a^  ■ 21.8  A.  The  potential  energy 
n(r)  takes  into  account  the  interaction  of  the  electron  with  the 
impurity  ion  and  with  the  image  of  the  impurity  ion  and  also  the 
screening  of  the  impurity  ion  potential  by  electrons  in  the  in- 
version layer.  In  the  present  paper,  we  neglect  screening  and 
take  U(r)  to  have  the  Coulomb  form 


U(r)  - - f 


where  r ■ (x^  + y^  + z^  . The  term  Involving  6 in  Eq.  (1)  repre- 


sents the  interaction  of  the  electron  with  its  own  image.  Our 


We  shall  be  particularly  interested  In  the  regime  of  large 
electric  fields  where  the  electric  sub-band  energies  are  large 
compared  to  the  impurity  binding  energy.  The  problem  is  then 

(7  'i 

analogous  to  the  problem  of  an  impurity  in  a high  magnetic  field^ 

and  the  problem  of  separating  the  electronic  and  nuclear  motions 

f 8 ^ 

in  the  theory  of  molecules  and  solids.  ' We  use  the  same  sort 
of  adiabatic  approach  that  has  proven  useful  in  the  latter  two 
problems . 

Let  us  write  the  Hamiltonian  in  the  form 

H - (3) 


where 


a? 


Ho  " “ Y fr  z 


z a 0 


H, 


U(r) 


(4a) 

C4b) 


Let  the  eigenfunctions  of  be  denoted  by  present 

paper  we  asstime  that  the  wave  function  i|r(r)  consists  of  a product 

_ 

♦ (r)  - X (x.y)f,  (z)  (®) 

^o 

where  f (z)  is  the  electric  sub-band  wave  function  of  lowest 
^o 

energy.  For  f (z)  we  choose  the  variational  form  of  Fang  and 
(9  ) 

Howard^  ® ' : 

f.  (z)  - (b^/2)^e“*’*/^z  , z > 0.  (6) 

For  the  lowest  sub-band  there  is  a series  of  bound  states. 

We  shall  discuss  the  ground  state  and  the  lowest  excited  states 
that  are  coupled  to  the  ground  state  by  electric  dipole  transitions. 


6. 


The  overall  wavefunctions  for  these  states  have  the  symmetries  of 


hydrogenic  2Pq  and  3d^^  states,  respectively.  For  these  states, 


we  choose  the  x (x,y)  to  have  the  variational  forms 


i -aop/Z 
- (a2/2rr)^e 


(Xo) 


4 i 

(a2/12TT)^(x  ±iy)e 


where  p 


+ y 


In  the  spirit  of  the  adiabatic  approximation,  one  would 


determine  the  parameter  b by  minimizing  the  expectation  value 


of  We  shall  generalize  the  situation  by  determining  all 


variational  parameters  through  minimizing  the  expectation  value 


of  the  total  Hamiltonian  H.  This  procedure  leads  to  different 


values  for  b for  the  2p^  and  3d  - . ...  . . ^ . 

o i;l  states  which  we  designate  by 


b2  and  b^,  respectively. 


III.  RESULTS 


A.  Energies 


The  evaluation  of  the  expectation  values,  £«  and  E_.  , 

3ci . 1 


is  straightforward.  The  values  of  the  variational  parameters  a^, 
in  each  case  were  obtained  from  the  equations 
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In  order  to  obtain  the  binding  energies,  we  must  obtain  the 
expectation  value  of  the  Hamiltonian  without  the  impurity  term. 

It  is  sufficient  to  consider  the  Hamiltonian  H^  given  by  Eq.  (4a) . 
We  use  the  variational  trial  function  given  by  Eq.  (6) • The 
expectation  value  of  the  Hamiltonian  H^  calculated  with  this 
trial  function  is  given  by 


" i yD  + i 6b  + 

S 


from  which  one  determines  the  optlmvun  value  of  b by  minimizing 

ik 

E . From  the  minimized  values  of  E-  , E.,.  , and  E_,  the  binding 

s 2Pq  s 

energies  are  calculated  from  the  equations 


■ ®2Po 


( 10a) 


®3d^l  - - =3d^i 


(10  b) 


Ve  tiave  calculated  the  variational  parameters  and  the  binding 
energies  for  the  2p^  and  3d^j^  states  for  electric  fields  ranging 
from  zero  up  to  lO^^esu.  The  results  are  tabulated  in  Table  I. 

The  binding  energies  are  plotted  as  functions  of  the  electric 
field  in  Fig.  1.  We  see  that  the  binding  energies  Increase 
fflonotonically  with  increasing  electric  field  and  approach 
asymptotic  limiting  values  at  very  large  fields.  The  limiting 
values  at  high  fields  correspond  to  the  binding  energies  of  two- 
dimensional  hydrogen-like  atoms.  The  binding  energies  of  the 
2p^  and  Sd^^^  states  increase  by  factors  of  over  ten  and  four, 
respectively,  from  zero  field  to  very  high  fields.  Associated 
with  the  Increases  in  binding  energies  are  decreases  in  the 
effective  Bohr  radii  in  the  plane  of  the  oxide-semiconductor 
interface.  The  effective  Bohr  radii  are  inversely  proportional 
to  the  variational  parameters  a.2  ^3  decrease  by  factors 

of  over  four  and  two  in  going  from  zero  field  to  very  high  fields 
for  the  2p^  and  3d^j^  states,  respectively. 

We  have  also  calculated  the  binding  energies  as  functions 
of  the  distance  of  the  impurity  from  the  silicon-oxide  interface 
when  the  impurity  is  in  the  oxide.  The  results  are  given  in 
Figs.  2 and  3 for  several  electric  fields.  We  see  that  the 
binding  energies  decrease  with  increasing  distance  of  the 
Impurity  from  the  Interface. 
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IV.  DISCUSSION 

In  our  previous  paper on  this  subject,  we  employed  the 
variational  trial  functions  of  Karpushin^^^  and  Bell  et  al^^^ 
to  investigate  the  effect  of  an  electric  field  normal  to  the 
Si02  ” Si  interface.  These  trial  functions  do  not  have  the 
form  specified  by  Eq.  (5)  and  do  not  give  the  correct  limiting 

4c 

values  for  the  binding  energies  at  high  fields  - namely,  4.0  Ry 

for  the  2p  state  and  0.444  Ry  for  the  3d^,  states.  The  trial 
o ±1 

functions  used  in  the  present  paper,  on  the  other  hand,  give 

the  correct  limiting  values  of  the  binding  energies. 

These  limiting  values  are  the  binding  energies  of  the 

appropriate  states  of  a two-dimensional  hydrogen-like  atom. 

It  should  be  noticed  that  the  two-dimensional  value. is  closely 

approached  for  the  2p^  state  only  at  very  high  electric  fields 
9 10 

in  the  10  -10  esu  range.  Such  fields  are  in  fact  not  realizable 


in  the  SiOg  - Si  system  since  the  energy  barrier  between  the 

SiOg  and  Si  is  only  -3  eV.  Physically  achievable  fields  may 

3 5 

reach  as  high  as  10  esu  (3  x 10  volts/cm)  where  the  2p^  state 
binding  energy  is  intermediate  between  the  two-dimensional  and 
three-dimensional  limits.  Thus,  the  use  of  the  two-dimensional 
binding  energy  in  discussing  the  properties  of  inversion  layers 
may  be  somewhat  misleading. 

The  calculations  reported  in  this  paper  are  based  on  a 
nvunber  of  simplifying  assumptions.  We  have  already  alluded  to 
our  use  of  an  infinite  potential  barrier  at  the  Si02  - Si  inter- 
face. It  would  be  desirable,  particularly  at  the  higher  electric 


Let  us  now  compare  our  theoretical  results  with  the  experi- 

( 6^ 

mental  results  of  Hartstein  and  Fowler'^  who  have  measured  the 
conductivity  of  n-channel  silicon  MOSFET  devices  in  which  Na^ 
ions  were  diffused  through  the  oxide  to  the  oxide-semiconductor 
interface.  Over  a certain  temperature  range,  the  peak  conduc- 
tivity can  be  described  by  an  activation  energy  which  Hartstein 
and  Fowler  identify  as  the  binding  energy  of  an  impurity  band. 

This  binding  energy  for  the  case  of  a narrow  impurity  band  should 
be  identifiable  with  the  binding  energy  of  an  isolated  impurity 
that  we  have  considered  in  the  present  paper.  By  varying  the  sub- 
strate bias  Hartstein  and  Fowler  were  able  to  vary  the  electric 
field  in  the  inversion  layer  and  thus  were  able  to  study  the 

activation  energy  as  a function  of  electric  field.  They  found 

11  -2 

that  in  a sample  with  5 x 10  cm  oxide  charge  density  the 
activation  energy  varied  from  18  meV  at  zero  substrate  bias  to 
25  meV  at  -15V  substrate  bias.  The  corresponding  electric  fields 
are  estimated  to  be  19.9  esu  and  60.4  esu,  respectively.  From 

Fig.  1,  we  see  that  our  calculated  binding  energies  for  the  2p^ 

* 

ground  state  at  these  two  electric  fields  are  0.50  and  0.60  Ry  , 
respectively,  or  21  and  25  meV,  respectively.  Thus,  our  theo- 
rectical  results  agree  rather  well  with  experiment. 

In  fact,  the  good  agreement  with  experiment  is  probably 
fortuitous  because  of  the  many  approximations  we  have  made. 

For  example,  the  impurity  ions  are  not  exactly  at  the  interface, 
and  the  Interface  Itself  is  not  a geometrical  plane  as  we  have 
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fields,  to  use  a finite  barrier  of  ~ 3 eV.  One  would  then  have 
to  allow  penetration  of  the  electron  into  the  oxide.  In  the 
present  work  we  have  assumed  that  the  impurity  ion  is  located 
exactly  at  the  interface.  In  an  experimental  situation,  however, 
the  impurity  will,  in  general,  be  situated  at  some  finite  distance 
from  the  interface.  This  possibility  is  taken  into  account  in  the 
Stern-Howard  formalism  and  is  currently  being  investigated  by 
detailed  calculations.  Another  point  which  should  be  raised  is 
that  the  interface  is  not  a mathematical  plane,  but  in  fact  is  a 
somewhat  diffuse  entity  extending  over  a region  of  5-lOA.^^^  ^ 

This  might  be  taken  into  account  by  using  a position-dependent 
dielectric  constant  which  varies  smoothly  from  the  value  for 
silicon  to  that  for  Si02  across  the  transition  region.  Another 
assximption  which  we  have  made  is  to  ignore  the  presence  of  the 
constant  energy  ellipsoids  whose  major  axes  are  parallel  to  the 
interface.  Taking  into  account  these  ellipsoids  would  lead  to 
additional,  higher  lying;  bound  states  and  some  modification  of  the 
bound  states  already  considered  due  to  intervalley  coupling.  We 
have  also  assumed  that  the  electric  field  due  to  the  inversion  layer 
is  a constant,  whereas,  in  reality,  this  is  not  so.  Nevertheless, 
the  work  of  Stem  and  Howard'  indicates  that  the  variation  of 
the  electric  field  over  the  dimension  of  a bound  state  is  rather 
small.  Finally,  we  note  that  our  use  of  the  effective  mass  approxi- 
mation may  lead  to  some  error,  particularly  at  the  higher  electric 
fields  where  the  bound  state  wave  functions  are  compressed  very 
close  to  the  interface. 


I 
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assumed,  ffe  have  neglected  whatever  screening  and  impurity 
banding  are  present  in  the  experimental  situation.  Also,  at 
the  relatively  low  electric  fields  in  the  experiments,  our 
variational  wave  function  can  probably  be  improved.  All  of 
these  aspects  are  now  being  investigated. 
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Figure  Captions 

Fig.  1.  Binding  energies  versus  electric  field  for  the  2p^ 
and  3d^^  states.  The  dashed  lines  are  the  two- 
dimensional  values . 

Fig.  2.  Binding  energy  of  the  2p^  state  versus  distance 
of  the  impurity  from  the  interface  for  several 
electric  fields. 

Fig.  3.  Binding  energy  of  the  3d^  states  versus  distance 
of  the  impurity  from  the  interface  for  several 
electric  fields. 
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